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Introduction 

Strictly  speaking,  a  potential  funculon  for  a  molecule  Is  purely 
a  theoretical  intermediate  in  the  approximate  calculation  of  Che 
energy  levels  and  not  an  experimental  observable.  There  would  seem 
to  be  little  reason,  then,  to  speak  of  an  "experimental”  potential 
function.  On  the  other  hand.  It  is  Impractical,  If  not  Impossible,  to 
compute  with  spectroscopic  accuracy  the  theoretical  potential  functions 
for  all  molecules  of  physical  Interest.  Thus  It  is  not  surprising  that 
numerous  methods  have  been  developed  for  estimating  the  potential 
function  from  experimental  data.  A  good  summary  of  these  methods  may 
be  found  In  the  review  article  by  Varshnl.^ 

As  no  two  of  these  methods  yield  the  same  potential  function,  It 
Is  of  some  interest  to  know  as  much  as  possible  about  the  theoretical 
potential  function  itself.  Again,  Varshnl  has  summarized  several 
basic  conditions  on  the  exact  potential  function.  There  are,  however, 
probably  a  wide  variety  of  theorems  which  can  be  proven  concerning 
theoretical  potential  functions.  Most  of  the  theorems  presented  in 
this  paper  take  the  form  of  Inequalities  which  place  rather  weak 
restrictions  on  the  potential  tunction. 

Discussion 

If  H^(R)  denotes  a  set  of  electronic  Hamiltonians  labeled  by 
the  parameter  R  and  E^(R)  denotes  the  corresponding  energy,  then  In 
atomic  units, 


-2 


H  -  T  +  V 
e  e  e 


where 


r  ^  ^  ^  r- 


The  potential  function  B(R)  is  then  given  by 


E(R)  -  E  (R)  +  z.z. /R  .  . 
e  a  b 


From  the  fact  that  T  and  V  are  homogeneous  functions  of  the 

e  e 

unit  of  length,  it  follows  that^T^/^R  -  -2T  /R  and -^V  /^R  ■  -V  /R  . 
This  yields  the  well  known  equations  , 


R^H  /■^R  -  -2T  -  V 
e  e  e 


-Tg  -  R  +  B  (RHg)/^  R  (2b) 

-  RdH  /3R  +  2H  -R"^^(R^H  )/3  R.  (2c) 

6  6  6  6 


Another  relationship  which  is  not  so  well  known,  can  be  found  from 


R^  /^R^  -  6T  +  2V 
e  e  e 


by  substituting  for  and  to  obtain 


R^9^H  +  4R9H  /3R  +  2H  -  0 

e  e  e 


(3) 
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If  the  symbol  )  E^R'>  is  used  to  denote  an  eigenstate  of  H^(R), 
then  Eg°(R)  “  *g°R^  where  1^°  Is  the  lowest  elgen> 

value  of  H  •  Since  H  (R)  has  no  singularities  in  R,  it  may  be 

G  G 

expanded  in  the  form 


H  (RJ  -  H  (R)  +  (R  -R)  ^  H  (R)/^  R  +  7(R,-R)^  V  (R)/^  R^+0(R, -R)^. 
eie  i  e  zi  e  i 


By  perturbation  theory  it  follows  at  once  that, 


Urn  B  “(RJ  “  8  °(R) 
R£^  R  ®  ® 


(5a) 


lim  dE  °(R  )/dR  3  dB  °(R)/dR  -  < E^°R  j  ^R  I  E  °R>  (5b) 

6li  6  G  G*G/^ 

1 

lim  d^E  °(R,)/dR.^5  d^E  °(R)/dR^-<B  °R  B  °R> -2R’^f (R) 

^  t>  Gil  G  G  G  G^ 


Rj-^R 


where  f (R) 


r2  O  I  V>/ 


E  >E 
e  e 


E  -  E 

e  e 


(5c) 


>  0.  (6a) 


The  formula  for  f(R) 


can  be  simplified  by  equation  (2a)  to 


Since  the  coordinate  transformations  which  commute  with  H  also 

e 

commute  with  and  separately,  it  follows  that  the  sum  over 

B^  ^  E^°  really  only  extends  over  states  of  the  same  symmetry  as  the 
ground  state. 

Equation  (5a)  is,  of  course,  the  Hellroan-Feynman  theorem 

dE  /dR  *  H  R^  .  Equation  (5b)  leads  to  the  conclusion  that 
2  ®  2  y  2  *  2\. 

d  E^/dR  K  H^/^S  /  for  the  ground  state  where  the  bracket 
denotes  the  average  value. 
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A  result  similar  to  this  can  be  obtained  for  an  optimized  para- 
meter  wave  function  ^  ^l’^2  *”  ^N’  *^1  where  the  a^  are 

variational  parameters  and  has  no  explicit  R  dependence.  If 
E  ~  ((p,  H<l>)  /  (p.p  ),  then  by  an  optimized  wave  function  it  is 
meant  that  "SE/^a^  ■  0  for  all  i  and  R.  But 

M 

di/dR  -  (S?/aR)^  ^  +y  (da,/dR)(^E/3a,)  (7a 


4  -  ^  [y 

^  dR  \  -^R/  dR  V  dR  ' 


Since 


i-  /ii).  ix  *  yi^)  ^  (^) 

dR  \^rJ  4-_\dR/  ■^a^S'SR/ 


it  follows  from  equation  (7a)  and  ■  0  that 


_i/M\ 

dR  \  ^  R/ 


Z  (■d^'da)  1^^ 


Equation  (7b)  then  becomes 


ii  _  a  ^E 


.  Y 

^  \  dR/l  dR/'^ajJaj 


But  the  condition  that  E  be  a  minimum  with  respect  to  all  the  a^ 
is  precisely  that  the  quadratic  form 


I  V  P. 

^  i  J 


« 


be  positive  definite.  Hence  it  follows  that 


Thus  for  either  an  optimized  approximate  solution  or  an  exact  solution 
there  is  an  inequality  of  the  form, 

d^E^/dR^<  .  (10) 

2  2 

If  the  expression  from  equation  (3)  is  substituted  for  ^  H  , 

2  2  2  ® 
this  inequality  becomes  d  (R  B  )/dR  <  0,  or  in  the  case  of  the 

2  2  2  ^ 

exact  energy  d  (R  E^)/dR  =  -2f(R),  Another  inequality  is  readily 
obtained  from  the  fact  that  is  positive.  Hence  the  inequalities 


2  9  2  2  2  2 

d  (R  E)/dR  S  d  (R  Eg)/dR^ 

< 

0 

(11a) 

d  (RE)/dR  c  d  (RE  )/dR 

e 

< 

0 

(11b) 

E  -  z  z^/R  =  E 

a  b  e 

< 

0 

(11c) 

must  be  satisfied  by  the  ground  state  potential  energy  of  every  dia¬ 
tomic  molecule. 

The  first  of  these  inequalities  places  a  restriction  on  the  force 

constant  of  the  form, 

*  • 

r^^  k  <  -  2  B(r)  (12) 

Since  the  zero  of  energy  is  fixed  by  equation  (1)  at  the  energy  of 

the  system  when  each  electron  and  nucleus  is  isolated  from  all  others, 

this  Inequallky  is  much  too  weak  to  be  useful.  It  is  interesting  to 

4 

note,  however,  that  Sutherland  has  found  a  closely  related  quantity, 

r  \  /D  ,  to  be  approximately  constant  within  the  same  class  of 
e  e  e 

molecules.  The  second  inequality  merely  restricts  the  slope  dE/dR 


-6 


since  it  implies  R  dE/dR  <  -  E  . 

The  Importance  of  the  inequalities  (11)  lies  in  their  implications 

concerning  the  smoothness  of  E(R).  Integration  of  these  inequalities 
2  2 

shows  that  RE,  R  E  >  and  d(R  E)/dR  must  all  be  monotonic  decreasing 

®  2 

functions  of  R.  Also,  RE  is  monotonlc  Increasing  to  a  positive 
maximum  and  monotonic  decreasing  beyond  this  maximum. 

Repeated  integration  of  inequality  (11a)  yields  the  inequality 


R.^E(E.)  -  R,^E(R,)  R,^E(R,)  -  r/e(R  ) 

V-R^,  ^  ^  '  l  -  'R  '  for  R,?R,>,R,?R^.  (13.) 

4  3 


^2  “  ^ 


This  should  be  a  useful  checking  formula  for  an  empirical  potential 
function  if  the  four  points  are  closely  spaced.  For  R2"R2"R  this 
may  be  written  as 


Pf 


\RiCRiE(Ri>  - 


h  -  ^ 


for  R,^  R^  R, 


(13b) 


Sab 

This  states  that  a  Fues  potential  ,  g2  “  ^  »  which  is  correct  for 

two  values  of  R  will  be  too  low  for  intermediate  values  of  R. 

2  2  2 

Double  integration  of  d  (R  E)/dR  ■  -2f(R)  yields 

Rj^E(Rj)  =>  R^E(R)  +  (Rj-R)  [^2RE(R)  +  R^dE(R)/d^  -  (Rj-R)^f  (^  ) 

for  Rj  ^  f  ^  R  (14) 

Some  properties  of  the  function  f(R)  can  be  obtained  from  the  limit- 
2  2  2 

ing  values  of  d  (R  E)/dR  .  For  instance, 

lim  d^(R^Ee)/dR^  -  2E  (*>«)  -  -2f(«*o). 

R-^**  ® 
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Thus,  £(oo)  is  the  negative  of  the  energy  of  the  separated  atoms. 
Also, 

lim  d^(R^E  )/dR^  -  2E  (0)  -  -2f(0). 

R-0  ®  ® 


So  £(0)  Is  the  absolute  value  of  the  energy  of  the  united  atom.  This 
gives  the  Interesting  'theorem, 

K-  ■  -  O  _  - = - STJ —  -  -  W  .  - y-  '  ,  fl —  (15) 


a  »  „  o  E  — E 

*  • 


,rv  o  *  -  E  ‘ 
E  >B  a  a 

a  a 


where  and  are  the  kinetic  and  potential  energy  operators  for 

an  atom  and  the  j  are  the  energy  eigenstates  of  that  atom.  It 

also  follows  from  this  that 


O  I  ..  2  ,  .  o. 


(16a) 


o  I  -  2  I  _  o 


(16b) 


which  gives  a  lower  bound  to  the  mean  square  deviation  of  and 

for  an  atom  in  its  ground  state. 

It  is  clear  from  this  that  £(R)  Is  rather  large.  In  fact,  for 
the  hydrogen  molecule  in  Hartree  atomic  units, 

f(0)  -  2.904  ,  f(rg)  -  0.808  ,  f(e«)  -  1  . 

For  other  molecules  f(R)  can  be  expected  to  be  even  larger.  Hence 

2 

the  last  term  in  equation  (14)  cannot  be  neglected  unler  (R^-R) 
is  small  compared  to  (R^-R). 

If  the  Inequality  (11a)  is  written  in  the  form 


-8- 


and  equation  (2b}  is  substituted,  the  inequality 

» 

d  T  /  dR  >  -  2  T  /R 
e  '  e 

is  obtained.  Or  in  a  different  form, 

d  Tg/dR  -  ^dTjdKy  =  -r‘^  d^(R^E)/dR^  >  0  (17a) 

d  V/dR  -  <dV/dR>  -  +r"^  d^(R^E)/dR^  <  0  (17b) 

Thus  it  is  seen  that  dE/dR  is  obtained  correctly  from  ^  dH/dR^ 
because  dl/dR  has  been  underestimated  by  the  same  amount  that  dV/dR 
has  been  overestimated. 

Since  the  derivation  of  the  inequalities  (11)  holds  equally 

well  for  the  lowest  excited  state  of  each  synmetry  it  follows  that 

the  inequalities  (12),  (13),  (14)  and  (17)  apply  to  the  potential  • 

o 

curves  for  these  states  also.  If,  in  equations  (15)  and  (16),  E 

1  * 

and  E^  are  Interpreted  to  mean  the  first  two  energy  levels  of  a  ^ 

given  symmetry,  equations  (15)  and  (16)  can  also  be  applied  to  excited 
states. 

Conclusion 

2  2  2 

The  inequality  d  (R  E)/dR  ^  0  places  certain  rather  weak 
restrictions  on  the  form  of  the  potential  curve  for  the  ground  state 
of  a  diatomic  molecule.  These  restrictions  probably  do  not  Justify 
the  amount  of  smoothness  usually  assumed  in  constructing  an  approxl* 
mate  potential  curve  by  empirical  methods  but  they  do  imply  that  the 
exact  potential  curve  is  not  too  badly  behaved. 


« 
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